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Abstract 

We construct lattice parafermions for the Z{N) chiral Potts model in terms of quasi-local 
currents of the underlying quantum group. We show that the conservation of the quantum 
group currents leads to twisted discrete-holomorphicity (DH) conditions for the parafermions. 
At the critical Fateev-Zamolodchikov point the parafermions are the usual ones, and the DH 
conditions coincide with those found previously by Rajabpour and Cardy. Away from the 
critical point, we show that our twisted DH conditions can be understood as deformed lattice 
current conservation conditions for an underlying perturbed conformal field theory in both 
the general N > 3 and N = 2 Ising cases. 


1 Introduction 

The chiral Potts model is a two-dimensional statistical model defined by spin variables subject 
to a Z^v-symmetric, local interaction, and was introduced in the 1980s as a lattice model for 
commensurate-incommensurate phase transitions [T]. A few years later, it became of great 
interest for mathematical physics as a solution of the star-triangle equations [2] in which the 
Boltzmann weights do not satisfy the difference property, and also as a superintegrable [3] 
generalisation of the Ising model. An important step in the understanding of the model was 
its identification [3] as a descendent of the well-studied six-vertex model - more precisely, the 
integrable chiral Potts model is based on Zjy cyclic representations mm of the affine quantum 
algebra Uj(s[ 2 ) (see also [7]). This class of representations was then studied in detail and 
generalised to other quantum algebras in [EHio]. 

The chiral Potts model displays very peculiar physical features. In the superintegrable 
regime, it describes a commensurate-incommensurate phase transition in an intrinsically aniso¬ 
tropic lattice model This behaviour is believed to appear also in the ordinary integrable 

regime, where the chiral Potts model provides an integrable chiral deformation m of the Fateev- 
Zamolodchikov (FZ) clock model [16]. The latter is an integrable ZAr-symmetric spin model 
whose scaling limit is described by the ZAr-parafermionic current algebra IIZ!, an extension of 
the Virasoro algebra generating the spectrum of a Conformal Field Theory (CFT). By a simple 

*ikhlefSlpthe.jussieu.fr 

^R.A.WestonShw.ac.uk 


1 



inspection of the Boltzmann weights, Rajabpour and Cardy have observed [TB] that one of 
the parafermionic cnrrents {ipk} of the Zjv-parafermionic CFT has a lattice analog in the FZ 
clock model, which is a discretely holomorphie operator, i.e., it satisfies a discrete version of the 
Cauchy-Riemann equations. Discretely holomorphie parafermions have been found empirically 
in a number of critical lattice models [MSI], and are a crucial ingredient to the rigorous study of 
bulk [22U26] and boundary m critical properties. In a recent paper [28], we have shown that the 
origin of discretely holomorphie parafermions for loop models can be traced to the underlying 
quantum algebraic structure, following the construction of quasi-local conserved currents by 
Bernard and Felder |29j . 

The object of the present paper is to use the integrability and Uq{5[2) symmetry of the 
chiral Potts model to (i) explain the algebraic origin of the discretely holomorphie parafermions 
observed [T8| for the FZ clock model, and (ii) extend the discrete Cauchy-Riemann equations to 
the chiral regime around the FZ point, and analyse their physical meaning in the scaling limit. 

The paper is organised as follows. Sections O [3] and H] are devoted to reviewing some back¬ 
ground material, respectively on the basics of the chiral Potts model, the Bernard-Felder con¬ 
struction, and the ?7q(sl2) symmetry of the chiral Potts model. In Section [5] we construct the 
quasi-local operators associated to generators of Uq{5{2) (which reduce to the lattice parafermions 
of [T8| at the FZ point), and we give an explicit form of the linear relations generalising the 
discrete Cauchy-Riemann equations for these operators. In Section [6] we interpret these linear 
relations in terms of perturbed CFT m and discuss both the general N >3 model and N = 2 
Ising case. Finally, we summarise our findings in Section [T] 

2 The Chiral Potts Model 

2.1 Definitions 

The chiral Potts (CP) model [TH3| is a statistical model where the variables are spins aj G Zjv 
living on the sites of a square lattice C. The Boltzmann weight of a spin configuration {aj} is 
invariant under a rotation of all spins {aj aj + 1 mod N), and is specified by the discrete 
functions VF{ij) associated to the edges of C: 


(ij) 


( 2 . 1 ) 


where (ij) denotes a pair of neighbouring sites, connected by an oriented edge i ^ j. Let us 
describe the specific choice of weight functions VF{ij) which renders the model integrable. In 
addition to the number of spin states N, we fix an external real parameter k' > 0. Each rapidity 
line carries a spectral parameter ^ given as a triplet of complex numbers ^ = (x, y, //) obeying 
the algebraic equations 

+ y^ = k{l + x^y^) and /i^(l - kx^) = k' , (2.2) 

where we have set k = y/\ — k'"^. A SW—)-NE (resp. NW—s-SE) edge crossed by rapidity lines 
(r, s) is assigned the weight function Wrs (resp. Wrs)^ definec|i| for a G {0,1,2,... } as 


Wrs{a) 



{lark's) X 

£=1 


XrOJ — XgU}^ 
ys - yrUJ^ 


(2.3) 


where u = exp{2iiT/N), and we have used = {xr,yr,kr) and = {xs,ys,ks) to denote the 
spectral parameters attached to the rapidity lines r and s, respectively. 

^ The conditions (EJI ensure that Wrs and Wrs are well-defined, i.e., Wrsia + N) = Wrs{a) and Wrs{a. + N) = 
Wrs{a). 
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The CP weights are represented by 


Wrsia-b) 




(2.4) 


We sometimes use an alternative graphical notation for CP weights that emphasises the fact 
that we can associate them with rhomboids on the covering lattice (i.e. the union of the CP 
lattice points, denoted by •, and the dual lattice points denoted by by o). This notation is 


Wrs{a — b) = am^ - , Wrs{a — b) 


A homogeneous chiral Potts model partition function Z^s is given as the sum over the height 
variables (i.e., the indices a E {0,1, • • • ,N — 1} at the positions marked by •) associated with 
a lattice with spectral parameters distributed as shown in Figure [H in which all diagonal 
lines have downward arrows, which we have omitted for clarity. 



Figure 1: A homogenous Chiral Potts Lattice 

The above weight functions satisfy [2] the star-triangle equations; 

N-l 

Wrsia - d)Wrtid - b)Wstid - c) = Prst X Wrs{c - b)Wrt{a - c)Wst{a - b) , (2.5) 

d=0 

for any fixed spins {a,b,c), and the overall factor prst is a function of ■ 



2.2 Crossing symmetry 

It is simple to see that the CP weights obey the crossing symmetry relations 


Wrs{a) 
where (x, y, p)* 


Ws*ria), Wrs{a) = Ws*r{-a) , 
{u~^y,x, 1/p), 


which can be indicated graphically by 



( 2 . 6 ) 
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2.3 Alternative parameterisation 


The spectral parameter ^ = (a:, y, y) subject to the conditions (|2.2I) can be conveniently repa- 
rameterised as 

(2.7) 


^ _ gi(M+ 0 )/V 




where the variables (0, 0, u) are now related by 


sin(/> = —/csintt, sin ^ cos tt, cos(p = k' cos( 

k' 


( 2 . 8 ) 


Note that (12.8p amounts to two independent relations, so that one complex parameter among 
((/>, (j), u) remains free. 

When approaching the self-dual line {(j) = /c' = 1}, it is convenient to scale u as 

TT 

u = -ilogk + — + u', 

where u' is finite. In the limit k' ^ 1, one then gets, on the self-dual line: 


sin (j) = sin cj) = -- 


2.4 Transfer matrix and spin-chain Hamiltonian 

In this paragraph, we consider a square lattice L tilted by 45°, so that the rapidity lines go along 
the horizontal and vertical directions, unlike in Figure [TJ We denote by Trs the transfer matrix 
comprising two horizontal rapidity lines with spectral parameter and 2L vertical rapidity 
lines with spectral parameter and we impose periodic boundary conditions in the horizontal 
direction. Then, as a consequence of the star-triangle equations (12.511 . the transfer matrices with 
different values of the horizontal parameter commute: 


[%i,s,Tr^A = ^ ■ ( 2 . 9 ) 

Moreover, when = Cs; the transfer matrix 'Trs reduces to a cyclic translation Hence, one 

can define the associated Hamiltonian in the usual way in the limit ^ by writing: 


Trs = e X {I - {ur - Us)'Hs + 0[{ur - n^)^]} . 


( 2 . 10 ) 


This results in the spin chain Hamiltonian acting on the tensor space Vl = ' 

times): 

L N-l 

«,, = T. —-EE ~ o '”1 


(L 


Otn — 


sin(7rn/A) 


z,r+ar,{x^x]^^r 

On = k' X 


N cos </>s ^ , 

J = 1 n=l 

exp [i{2n - N)(l)s/N] _ _ ^, exp [i(2n - N)^s/N] 


( 2 . 11 ) 


sin(7rn/A) 


with periodic boundary conditions = Xi. The operators Xj and Zj are defined through 

the elementary N x N matrices (A, Z) with coefficients 


^ab — ^ ^ab j 


Zab 


r(mod N) 
^a,b-l ’ 


( 2 . 12 ) 


which are characterised (up to a change of bases) by the relations X^ = Z^ = I and ZX = 
ojXZ. The operator Xj (resp. Zj) then acts as the matrix X (resp. Z) on the j-th factor of lA, 
and as the identity matrix on the other factors. Note that (12.lip is Hermitian for any choice of 
real parameters {(j)s,^s,k'). 
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2.5 Ztv charges and Kramers-Wannier duality 

Let R be the global rotation of spins: R = 0^=1 ^j- Since R and Rg commute, we can diagonalise 
them simultaneously. We denote the eigenvalues of R as w”™ with m ^ Zn, and we refer to m as 
the Ztv charge. We may also consider the Hamiltonian (j2.1ip with “twisted” periodic boundary 
conditions = uj^Xi, with m G Z^r, and we call fn the dual Z^r charge. 

In the present context, Kramers-Wannier (KW) duality amounts to the following non-local 
change of bases. We introduce for j G {1,..., L}, the dual operators 


z,+i/2 = y xj+i. x,+,/2 = nzt, 

£=1 

and we set by convention = I- We see readily that 

L N-l 

j=l n=l 


(2.13) 


(2.14) 


with = -^^+ 1/2 = ^5 Zjj_ii 2 X£j_ii 2 = X£j_ii 2 Z Hence we recover the 

original form (|2.11|) . except that the roles of an and have been exchanged. In terms of 
external parameters, KW duality acts as 


—> {(j),(l),l/k'). (2.15) 

Note that this transformation preserves the integrability condition k'coscj) = coscj). Let us 
examine its effects on the Z^r charges {m,rn). From the identities 

L L 

u^I=Ylzj = x\/2Xl+,/2 , I = = H , (2.16) 

i=i j=i 

it is clear that KW duality exchanges the two charges: 


(m, m) 


(m, m). 


(2.17) 


2.6 The Fateev-Zamolodchikov case 

When (p = (j) = 0 and k' = 1, the chiral Potts model reduces to the Fateev-Zamolodchikov 
(FZ) clock model |16) . In the scaling limit, the model is isotropic, and is described by the 
Z^v-parafermionic CFT [17] . 

The weights of the FZ clock model enjoy the difference property, i.e., Wrsia) and Wrs{a) 
are functions of Ug — Ur and a. Under these conditions, the star-triangle equations (|2.5p are 
consistent with the following embedding of the model in the complex plane: 


Wrs=*^ 





(2.18) 


where the angle 9 is defined as 

9 = Ug — Ur , (2-19) 

and the rhombi have a unit side length. Note that this choice also respects crossing symmetry, 
in that the embedding angle of the left hand side of the first crossing relation (12.61) is Ug — Ur, 
and the embedding angle of the right-hand-side is n — Ug + Ur (and similarly for the second 
crossing relation). 
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2.7 The Ising case 

The chiral Potts model with N = 2 corresponds to an Ising model with a partition function of 
the form 

{ij) 

where = Ki (resp. = K 2 ) if (ij) is a horizontal (resp. vertical) edge. In this case, 

the algebraic relations (12.21) can be parametrised by Jacobi elliptic functions of modulus km-- 


= —Vk sn/3 , y = -Vk- 


cn/3 


k- = 


Vy 


dn/3 ’ ^ dn/3 

For any value of k', the couplings are functions of j3s — /3r- 

= k'sc(i{K — /3s + Pr), = k'scd{Ps — Pr) 


( 2 . 21 ) 


( 2 . 22 ) 


where scd(u) = sn^/( cn^ dn^), and K is the complete elliptic integral of the first kind of 
modulus k. Using elliptic identities, one gets the relation 


sinh2ii'i sinh 2 iF 2 = 


1 

y 


(2.23) 


Expressions (12.221) show that the Ising model enjoys the difference property for any value of 
k. Requiring that the star-triangle relation can be represented as a geometric relation in the 
complex plane leads us to define the embedding angle (see Section ESI) as 


TT 


- /3r) • 


From (12.7p and (|2.21|) . we have the relations 


^ ik sn/3 cn/3 ^ i sn/3 dn/3 ^ ik' sn/3 


dn/3 


cn/3 


cn/3 dn/3 


(2.24) 


(2.25) 


In the critical limit fc —)• 0, a way to recover the expression (I2.19P for the embedding angle is to 
set 

K / \ 

(2.26) 


/3 = ^ (|logp + 2/3') , 


where /3' is finite \p is the nome of elliptic functions in (12.211) . with k ~ In this regime, 

we have the expansions of elliptic functions (see m)-- 


H{/3) = -iy^' + + 0(p3/2) ^ 

Hi{I3) = + 0(p3/2) ^ 

0(/3) = 1 - + 0(p3/2) ^ 

01 (/3) = 1 + + 0(p3/2). 


(2.27) 


Using these, we find that e®“ = —-|- 0(p^y), and hence ()2.24p reduces to (|2.19p up to 
corrections of order k^. 


2.8 Critical behaviour 

Returning to the general N > 3 case, we shall describe the critical behaviour of the chiral Potts 
Hamiltonian ()2.1ip in the space of parameters ((/), (^, A:®). 
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On each of the self-dual (SD) lines 

SDi = {(/) = (^, k' = 1} and SDa = {(^ = -0, k' = 1}, (2.28) 

the model is invariant (modulo a global spin reversal) under KW duality, and thus it is massless. 
On each of the SD lines, for cp ^ 0, the model is in an “incommensurate state”, (i.e. the ground 
state is in a sector with non-zero momentum), and remains critical |T5], with the same critical 
exponents as in the isotropic Zjv-parafermionic OFT, but its correlations become anisotropic 
in Minkowski space. These SD lines correspond respectively to a perturbation of the Zjy- 
parafermionic CFT by an operator of conformal spin -|-1 and —1. 

As shown in mHH], in the plane </> = </>, the model remains massless in a limited region 
around SDi, and then it undergoes a commensurate-incommensurate transition to a massive 
phase. It seems plausible that these results also hold outside this plane, i.e., there is a massless 
incommensurate phase around the plane {k' = 1}, and the model becomes massive for \k' — 1| 
large enough. At the FZ point, the massless phase is “pinched”, and any small perturbation 
outside the plane k' = 1 develops a finite mass. 

3 Non-local Operators and Quantum Groups 

In this section, we review the picture of non-local operators arising from quantum groups that 
was developed by Bernard and Felder in [29]. In [28], these currents were used in a direct way 
to construct discretely holomorphic operators in dense and dilute loop models. In the case of 
chiral Potts, we will show in Section [5] that the currents constructed using the method of [29] 
split naturally into two half-currents which obey a discrete holomorphicity condition. 

The starting point of Bernard and Felder [29] is to consider a quasi-triangular Hopf algebra 
A (also known as a quantum group) defined in terms of a set of generators {Ta, ©“ft}) 
a,b = 1,2,... ,n that have the relations 

ea’’e% = Sa,c and e\eb^ = Sa,c (3.1) 

and the coproduct structure 

A{Ja) = Ja^l + ea^^Jb, A(0„^) = 0/ ® 0,^ A{Q\) = 0% ^ 0% , 

where A : A ^ A (8* A denotes the coproduct (for a gentle introduction to quantum groups see 
for example |32|). Note that we use the convention that repeated indices are summed over. The 
antipode and counit that complete the Hopf algebra structure can be found in [2^ . 

It is helpful to introduce a graphical notation for representations of A. We indicate a 
representation of A by a line and the action of the above generators of A on this representation 
by 



Adopting the convention that composition AoB means that B is above A, the inversion relations 
(EH) are then represented as 




A- 


and 


(3.2) 
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The action of the generators on tensor products of representations is indicated by 


A(Ja) 





Ja ® 1 ® Jfe 


A(0a') 



<8) ej’ (g) 


Denoting the above coproduct by A^^) : ^ ^ (g) ^ , it is then possible to define a coproduct 

A^^^ (with L terms in the tensor product on the right) recursively by 

^(m+i) _ (/\(g)I(g)I(g).I)A(”^). Then it follows that the representation of A^^\Ja) on an 

L-fold tensor product is indicated graphically by 




2=1 



The above graphics becomes more useful when combined with the standard graphics for the 
R-matrix. The R-matrix of a quantum group, i? : Ri (g) V 2 ^ V 2 <8> Ri, is a map between tensor 
products of representations Vi and V 2 that commutes with the action of the quantum group A. 
That is, we have RA{x) = A(x) R. We can represent the R-matrix graphically by 



where the arrows serve to orient the picture and we view the R-matrix above as acting from top to 
bottom. We shall generally suppress these arrows. The commutation relations RA{x) = A{x) R 
then have the following simple graphical realisations when x = Ja, and 0“^: 



We now wish to define a non-local operator ja{x,y) associated with the insertion of the 
operator Ja at a point (x, y) in a 2D lattice model and an attached tail made up of tensor 




products of QJ' and 0 “b along some path leading to a marked point on the boundary of the 
lattice. In order to do this, we again follow the approach of Bernard and Felder. Suppose we have 
a 2D lattice A consisting of 4-vertices at points p G Let A' denote the lattice consisting of 
the points r which are the midpoints of the edges of A. Then if V{r) denotes a A representation 
associated with the midpoint r, the R-matrix associated with the vertex p ^ A will be a map 


R{p) : V{fi) ® Vin) ^ V{r2) ® Virs) 



(3.5) 


where the r, G A' are the indicated four midpoints surrounding the point p (some of the V)?. will 
need to be isomorphic for the R-matrix to exist, which we assume to be the case). We can then 
define a vector space 

Fa = (g) F(r-'), 

rsA' 

and a linear operator B : Va ^ V\, hy 

B = ^R{p). 

peA 

The partition function of the vertex model with Boltzmann weights specified by the R-matrices 
is given by 

Z = Trv^{B), (3.6) 

and the expectation value of any linear operator O : Va ^ Fv is given by 

{0) = ^TvvAOB). (3.7) 

The above expressions (j3.6l) and (j3.7p may at first seem unusual - the partition function and 
correlation functions are more commonly given as the trace over a ID tensor product space that 
would correspond to the Hilbert space in the quantum statistical mechanics interpretation of 
the partition function. However, it is a useful, and simple, exercise to check that the partition 
function does always reduce to the standard ID trace. In contrast, the expectation value (O) 
may be written as a ID trace only in the special case when the operator O acts trivially at all 
r in Va except those along some ID line in the lattice. 

In this paper, we are interested in operators O that act non-trivially on a set of midpoints r 
that will wind along a path 7 from a marked interior point to a marked point on the boundary of 
the lattice. In this case the more general expression (|3.7I) will be required. To be more precise, 
we consider the operator j2if^ ■ Va ^ Va constructed by the insertion of the representation 
of Ja on V{r) at the midpoint r, and the insertion of a “tail operator” constructed from the 
insertion of ©a^ or 0 “^ along a line of midpoints specified by the path 7 that terminates at some 
fixed, but arbitrary, point on the boundary of the lattice. An example is shown in Figure [3l 

The commutation relations with the R-matrix expressed by (|3.3D and (13.411 have two imme¬ 
diate consequences for expectation values (jd(r)). The second relation (|3.4I1 implies that the 
expectation value is independent of the path 7 and will depend only upon the insertion point 
r and the fixed boundary point. Thus we will from now on drop the 7 path superscript on the 
current. The other commutation relation (j3.3l) implies that when inserted into an expectation 
value we have 

ja{n) - ja{f 2 ) - jain) + jaiu) = 0 , (3.8) 

where r) are the four edge midpoints points surrounding any vertex - as indicated in (13.5|) . After 
embedding the lattice into the complex plane it is this relation (13.8p that was interpreted as a 
discrete holomorphicity relation in several examples in the paper |28| . 
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Figure 2; The insertion points and path of a non-local operator 


4 Chiral Potts Weights and Representation Theory 

In this section we review the construction of chiral Potts weights in terms of the representation 
theory of Uq{5[2) [SHTO]. The notation we use is that of [9]. 


4.1 The quantum afRne algebra Uq{s\ 2 ) 

We begin by defining [/^(sb), which is an algebra over C generated by Ci, , Zi {i = 0,1), 

where satisfy the standard relation of the quantum affine algebra 17q(sl2), and zq and 

zi are two new central elements. The comultiplication of is chosen as 

A(ei) = ei®I + Ziti®ei, A{fi) = fi 0 t~^ + z:r^®fi, 

A{ti) = A{zi) = Zi®Zi. 


The representations relevant to the chiral Potts occur when q = — exp(i7r/iV) where N G 
{2,3,4, • • • }, and we shall fix q to take this value from now on. We also define oj = (f‘. These 
representations are Wdimensional cyclic representations denoted V^r' and parametrised by a 
pair of points (r, r') ^ Ck x Ck- Here Ck is the algebraic curve given by (x, y, z) G such 
that 


+ y^ = k{l + x^y^), 


= 


k' 


1 — kx^ 


1 — ky^ 
k' 


where k'^ + k''^ = 1. 

If r = (x, y, z) G Ck and r' = (x', y', z') G Ck then the representation V^r’ is given by 

T^rr'{ei) = - y')X , 'Krr'ih) = - x'y,y !), 

xx'y/i' 

T^rr'ieo) = - x') , 7rrr^{fo)= (^Z-Z\x, 

[q - 1) V ^ Coy J 

^rr' {to} 7 '^rr' i^o} ^0 • 

Here, the objects X and Z are N x N matrices, such that ZX = loXZ and X^ = Z^ 
this paper, we shall fix X and Z as 


/I 

0 

0 • 

• 0 

0 \ 


/O 

1 

0 ••• 

0 

0\ 

0 

UJ 

0 • 

• 0 

0 


0 

0 

1 0 


0 

0 

0 

a;2 . 

• 0 

0 

, Z = 












0 

0 

0 ••• 

0 

1 

Vo 

0 

0 • 

• 0 



VI 

0 

0 ••• 

0 

0/ 


(4.1) 


I. In 
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The constanlH cq also appearing in the above satisfies Cq = q^xx'/{yy'). 

4.2 The ^-matrix 

Consider now the R-matrix R{rr',ss') : Vrr' ® Vss’ Vgs’ ® Vrr' that obeys R{rr', ss')A{x) = 
A{x)R{rr', ss') for x G The approach of [9] starts with the Ansatz that this R-matrix 

is of the factorised form 


R{rr', ss') = Sr'siTr's' ® Trs)Srs' , 

where 

^ ® ? ^Ts • ^sv ^ ^rs * 

Then the relation 

R(r/, ss^)[7rr^/(g) 7rss/(A(x))] = [tTs^/ 0 7rrr'(A(x))]R(r/, ss^) 

is ensured if S and T satisfy the stronger ‘sufficiency conditions’: 

(g) IT ggt (A(x))] = [ 7Vgfy>f (S' ^sr (A(x))]5'„/, 

{Tr's' <g) 1)[ (S ST (A(x))] = K fg! (S ^sr (A(x))](Tr/5/ (g) 1), 

(1 (g) rrs)[7rr/s/(g) 7rsr(A(x))] = [vr^/s'® 7rrs(A(x))](l ® T^s), 

*S*7*/5 (g) (g) ^A^x))]. 

These sufficiency conditions were in turn found to be satisfied by the choice 

N-l 


(4.2) 


(4.3) 

(4.4) 

(4.5) 

(4.6) 


Srs{Vei ®Ve2) = Wrs{.ei - £2){Ve2 ® '^ei) , 


TrsVs — ^ ^ iTfs(n)u£—a ; 
a=0 


with the coefficients given by 

Vs X'pUJ 


1 ) iJjg yrp XgOJ^ 


Wrs{^) XrUJ — XsCj'^' 


Wrsin - 1) 


ys - yrOJ" 


where we have now switched to the notation a = {xa,ya,fJ'a) for a point a G C^- These are the 
precisely the chiral Potts model Boltzmann weights (|2.3D which can be found in [2]. Defining 
components R{rr', ss')“^ of the R-matrix by R{rr', ss'){va ®vi,) = J2 Ri^r', ss')“^(xd ® Vc), then 

c,d 

leads via (SSI) to the factorised expression 


R(rr', ss')cd = Wr's{d - c)Wr>s'{a - d)Wrs{b - c)iT„/(a - b ). 


(4.7) 


It is useful to introduce a version of the standard 4-vertex graphical notation for R-matrices 
that is modified to deal with the representation Vrr'- We indicate the identity acting on the 




representation Vrr' by a directed double line 


R{rr', ss') = 


, and the R-matrix R(rr', ss') by 



^The other constants appearing in Section 4 of [5] are here fixed as ci = I/cq and Ki = l/(g^ — 1) 
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The components R{rr', are then indicated by 



with the graphical conventions of Section [21 

5 Construction of non-local operators in the Chiral Potts model 


In this section, we consider the non-local operators ja{r) discussed in Section [3] associated with 
the cyclic representations of C/g(sl 2 ) introduced in Section jH 


5.1 The current jgo 

Let us first consider the current associated with the generator bq := to/o- We have the co¬ 
products 


A(eo) = eo (8) I + to-^o ^ ® ^0 , and A(to2Q ^) = to^o ^ ® ^o-^o ^ > (5-1) 

and thus bq and are respectively generators of type Ja and 0^ in the notation of Section [3l 
It follows from (|4.1I) that the action of these generators on the representation 14-r' is given by 

T^rr'ieo) = X [x~^ - y~^TTrr'{toZQ^)] , TTrr'= fr fr'Z~^ , (5.2) 

where fr '■= ——— • 

—qXrlJ^r 

We now wish to follow the approach of Section O and consider the non-local operator BQ{x,t) 
associated with the insertion of the appropriate representation of the non-local operator 


• • • (g) toZQ ^ (g) toZQ ^ (g) eo . 


(5.3) 


The position (x, t) will correspond to a CP site (x, t) (indicated by a • in Figured]) and we define 
Bo{x,t) such that the cq ia i|5.3l) acts on the representation associated with a pair of diagonal 
lines either side of the point {x,t). To make this definition more precise it is useful to modify 
the graphical notation of Section [3l To this end we introduce the following representation of the 
diagonal action of X and non-diagonal action of TTrr'itoZQ^) on the representation Vrr- 


r 


X ~ 



7r(j.j,/) (tQZQ 





t t 


It follows from (15.ip and (|5.2p that the action of eo(x,t) splits into two ‘half-currents’ which can 
be represented graphically as 


Bo{x,t) = 



(5.4) 


Here (x, t) is the CP site marked by a ■ and the left tail will wind through the rest of the CP 
lattice. Up to considerations of boundary conditions, which we will not need in this paper, the 
position of this tail is arbitrary due to the commutation of t^ZQ^ (g) with the R-matrix. 
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Let us now consider the sufficiency condition (14.31) in the case when x = cq, r' = r and s' = s. 
Equation (|4.3p becomes 

TTfr (S' '^ss (A(e-o))] = [ TTgr ^ '^sr iA{eo))]Srs, 
which can be represented graphically as 


r s 



Several points are worth noting about this commutation relationship: 

1. The tail operator is naturally associated with the edges of the dual CP lattice with vertices 
indicated by o. 

2. We have appended a left horizontal tail to indicate that this relation may be embedded in 
the larger CP lattice with such a left tail. 

3. There is cancellation of four of the terms. 


After cancellation, we arrive at the four-term relation 



(5.5) 


Note that we have used the relation XZ~^ = q^Z~^X to rewrite the first term on the right- 
hand-side of this equation. 

In order to rewrite the relationship (|5.5p in terms of CP weights, we need to understand the 
effect of the tail operator purely in terms of the modification of CP weights. Recall that the 
action of the tail operator toz^^ on the representation Vrr' is given by 

TTrr'itoZQ^) = frfr'Z~^ , ( 5 - 6 ) 


where Z ^ \s the cyclic shift matrix that acts on canonical basis vectors Uq (a = 0,1, • • • , — 1) 

as Z~^Va = Va+i\ ^ • Hence, we can identify 


= ^Wrs{a-b + l), 


= ^Wrsia - b - 1), 


'\z\^ 


— frfs^^rsi^^ ^ “h 1)^ ^ 


(5.7) 


b 
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In this way we identify the tail operator as a disorder operator for the CP model. 

The diagrammatic identity of (15.5p can be written as a relation around any horizontal CP 
plaquette inserted into a larger partition function as 



Consider the following non-local operator jg^ (r) defined in terms of the above half-currents by 

jeo 2 > 


where 

- (T(ro-) corresponds to the insertion of X = ■ at embedded CP site rv 

- fieoir^) is the above tail/disorder operator ending at embedded dual CP site 

- T is ‘tail ordering’ defined as: 


quasi-local op. fj.eQ{f^)a{f„) with tail 


locally above fv if < Im(zo-) 

locally below fv if > Im(zo-) 


With this definition, the graphical relation (15.8p can be written simply as 

- Vr^jeoin) + 9^2/7^60(^ 2 ) - x-^jeoin) xj^jgoiu) = 0 , (5.9) 

where we have denoted the mid-edges of the plaquette as follows: 

ri^' ^\r4 

^ . (5.10) 


5.2 The operator and twisted Cauchy-Riemann equation 

Let us parameterise {x,y,fi) in terms of (n, (/>, i^) as in (12.7p . and introduce 

«! = -^-TT, 02 =- - -hVT, 03 = - - - , 04 =- - - . (5.11) 

The linear relation (j5.9p reads: 

e*('^’-+“^)/'^Je-o(n) - jg,{n) = 0 . (5.12) 


We choose 6 = Us — Ur ss the embedding angle in (12.1811 . so that the quantity aj coincides with 
the principal argument of (z^ — z^) on the corresponding edge of the plaquette. Then the linear 
relation takes the form of a “twisted Cauchy-Riemann” relation: 


i'Pr 


e N §zi OgQ^(ri) + eN dz 2 Cgo^(T 2 ) + e 7 5z3 + e 


is) I 
eo 


eo 


■^fc4Oe?(r4) = 0, (5.13) 


where we have introduced the lattice current 

Ogjif^ = exp[-iso(f)] jeo(r), (5.14) 

and where o(r) is the principal argument of the oriented edge (z^ — z^) carrying the point r, 
and the spin s is set to s = 1 — 1/iV. 
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5.3 The Wrs plaquette 

Equation (I5.13P is a discrete integral relation around a Wrs plaquette. There is a direct extension 
of the above arguments that leads to a discrete integral relation for Oe^iz) around a Wrs 
plaquette. The starting point is to consider the commutation relation (14.5p . that is 

(Trs (g) l)[7rsr <8* 7rsr(A(eo))] = [vr^s <8) vrsr(A(eo))](rrs ® 1) • (5.15) 

Using the splitting into half-currents given by dsap, this can be represented by 



Defining the operator (r) exactly as above, this leads to the discrete integral condition 


e "n SziO'g^ifi) + e 5z20g'J{f 2 ) + eN §z 30 gl\r 3 ) + eN = 0 , (5.18) 


around the embedded plaquette 


ri^ 

sC 

r2^ 


\r4 

\ 

> 

*r3 


(5.19) 


The coefficients appearing in (15.131) and (|5.18p are such that when we consider any larger region, 
the contribution from internal points cancel and we are left with a discretely holomorphicity 
relation expressed solely on the boundary of the region. 


5.4 Other currents 

It is possible to define quasi-local operators in terms of the half-currents associated with the 
other generators. We consider these in turn. 


5.4.1 The current for cq 

The coproduct and representation of sq are given by 

A(eo) = eo (8)I + io^o ® eo, A(to^o) = ® , (5.20) 

7r^^/(co) /?A Vr'^rr^ 5 ^rr^(^O'^o) ^ ? 
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where /3 := —q/{q^ 


1)^. Proceeding as above, may now define a new half-current 


Jeo((fV + f>)/2) = r(/reo(f^)cjt(f„)), 


where 


- (T^(fv) corresponds to the insertion of X ^ at embedded CP site fv 

- is a new disorder operator corresponding to the insertion of the operator t^ZQ along 

a path ending at the dual CP site r^. The effect of this disorder operator is similar to that 
of (j5.7p except that fr,s right-hand side. 

Defining = exp[-t-fsa(f)]jeo(f^)) we can follow the previous analysis to arrive at the 

discrete antiholomorphicity conditions 


e + e Sz 20 ^^^^(r 2 ) + 

e“^ dziO^~^^(fi) + e“^ dz 20 ^~^^(f 2 ) + e~ 


e^Sz 30 t\r 3 ) + e^ 
^6z30tHr3) + e-^ 


<5z4O(;*)(f4) = 0, 
Sz 4 O(;^^(n) = 0 , 


around the W plaquette (I5.10D and W plaquette (15.191) respectively. 


5.4.2 The current for ei 

A similar consideration leads us to define 

Jgi((fV + f>)/2) = T(fisi(r^)cT^f^)), 


where 

- (T^(ro-) corresponds to the insertion of X~^ at embedded CP site Po- 

- is a new disorder operator corresponding to the insertion of the operator 
along a path ending at the dual CP site r^. Noting that 

TTj-r' (tlZ^ ) — f f/ ^ 

Jr Jr 

and comparing to equation (|5.6I) . we see that this disorder operator has the same 

action as in (1521), but now with the arrow is directed outwards from the dual CP site r^. 

Defining Og^\f) = exp[—fsa(r)]jg^ (r), we obtain the W and W discrete holomorphicity condi¬ 
tions 






' 'n 5ziOgl\fi) + e "n 5 z 20 gi\f 2 ) + 6 ^ 5 z 30 gl\r 3 ) + 6 ^ 5 z 40 gl\u) = 0, 


e "n + eN 6 z 20 gl\f 2 ) + e ^ dzsO'f^’ifs) + e ^ [r^) = 0 . 








(5.21) 

(5.22) 


5.4.3 The current for ei 

Finally, by considering the half-currents that make up ei, we arrive at the following definition 
of a quasi-local operator 

jei((P<7 + P/.)/2) = r(^ei(P/.)cT(f<^)) , 

where 

- C7{f(j) corresponds to the insertion of X at embedded CP site Po¬ 


le 







- /Xei(^/i) is a disorder operator corresponding to the insertion of the operator tizi along a 
path ending at the dual CP site r^. This effect of the disorder operator is given by (|5.7I1 
with fr^s on the right-hand side and with the arrow leaving f^. 

Defining Oe^^\r) = exp[isQ;(r)]jei we obtain the W and W discrete anti-holomorphicity 
conditions 


e ^ + e~~^ S220e7^H'^2) + 


^fc-30(r^(f3) + e-^ 
e^ 6 z 30 (-^HrV + e^ 


Sz4O(;^^(f4) = 0, 

Sz4O(;^Hr4) = 0. 


In summary, Cj yield half-currents with discrete holomorphicity and spin s = (1 — 1/A^), 
whilst Ci yield ones with discrete antiholomorphicity and spin —s. The coefficients in the discrete 
relations around either (Pl, 7 *2, rs, 74 ) plaquette are just ^^±i(f)s/N^ 

cases. 


6 Physical interpretation 

In this section, we discuss the physical meaning of the linear relations derived above for the 
operators Og/ and Oe^ ■ For simplicity, we restrict the discussion to the case of 0\7 and the 
linear relation (j5.13l) around a Wrs plaquette, but very similar results hold for the other cases. 
Also, in order to lighten notation, we drop the indices cq in jeo(r) or Ogg\f), and simply write 
j(F) and C>(^)(r). 


6.1 Discrete linear relation in the vicinity of the FZ point 

At the FZ point (7r,s = 4>r,s = 0, ()5.13l) takes the simple form 

SziO^^^fi) + 6z20^^Hf2) + SzsO^^^fs) + 6z40^^\f4) = 0. 


( 6 . 1 ) 


Note that at this point 0^^\z) coincides with the lattice parafermion of |18] . and (16.ip is the 
discrete Cauchy-Riemann relation of the form d'ljjs = 0 found empirically in |18] . 

In the vicinity of the FZ point, using (j2.8l) . we can write the linearised relations 


(/)r = cos 6 (ps + i sin 9 (ps + 0{(p^, (p ^), cpr = i sin 9 (pg + cos 9 (ps + 0{(p^, (p^). (6.2) 

Hence, if we introduce the notation (p^ = {(pzt 4>)/2, we get (pf- ~ e^^'^cpf, and (I5.13p becomes 
6ziO^^\ri) + 6z20^^\r2) + SzaO^^^fs) + Sz^O^^^n) = 


— a 


+ t-^o^^\f2) + + r^oW(r4) 

+ a7 (fi) + (f2) + (r3) + (f4) 

+ i{atf + t-^0^^-^\f2) + 

- i{a;f + tO^^-^\f2) + 

+ 0{{aff), 


where we have set 


t = —ie*® 


N 


e 


N 


(6.3) 


(6.4) 
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Note that the sum of coefficients multiplying the O's in each bracket on the RHS of (j6.3p is 
4sin0, which is proportional to the area of the plaquette. In the case 0 = tt/2, then t = 1 and 
takes the simpler form: 

6 ziO^^\ri) + 5z20^^\r2) + SzsO^^^ifa) + fc40^*^(r4) = 

- + a-d(^-^)(p) + i[(a+)2 - + 0((a±)=*), (6.5) 

where we have defined p as the center of the plaquette, and 

d(p) = o(fi) + o(f2) + o(f3) + oiu) ■ 

6.2 Comparison with perturbed CFT 

In general, consider the perturbed action [30] 

S = Soft + g j (fr > 

where Scft is the action of some CFT, and ^ is an operator in the spectrum of this CFT. 
Suppose the unperturbed theory possesses a holomorphic current with conformal spin s, 

and moreover, assume an OPE between the current and the perturbing field of the form: 

= ••• + +■■■ ( 6 - 6 ) 

By dimensional analysis, x must have conformal weights = s + h — 1 and hy. = h. In 
the perturbed theory, 'ipg is no longer a conserved current. More precisely, using the identity 
dzjl/(zi — Z 2 )] = vr5(ri — r 2 ), one gets at first order in g: 

B'lpsiz.z) = Tig x{z,z). (6.7) 

So this simple line of arguments tells us how the conservation equation for the current ips is 
modified by the perturbation. 

Let us now specialise to the Zjv parafermionic CFT for N > 3 (the case = 2 is treated 
separately in the next section). The chiral Potts model corresponds [151133] to a perturbation of 
this CFT by the energy operator e, together with the leading spin ±1 operators: 

S = Sfz + j (fr[5+^+{z,z) + 5-^-{z,z) + Te{z,z)] . (6.8) 

The energy operator e has dimensions = he = 2/{N + 2), and the spin ±1 operators (which 
are actually the descendants W-ie and IT_ie in terms of the underlying Wn algebra) have 
dimensions (/i$_^,/i$^) = {he + 1,/ie) and (/i$_,/i$_) = {he,he + 1). The parafermion of Z^r 
charges m = m = 1 has conformal spin s = 1 — 1/A^, and its conservation equation is modified 
by the perturbation according to dEZ]): 

dil;s{z, z) = 116 + x+{z, z) + vrd. x-{z, z) + vrr xo(^;, z ), (6.9) 

where the operators on the right-hand side have Z^v charges m = m = 1, and conformal 
dimensions: 

{^X+ 1 ^X+ ) ~ {^£ S,he) , {f^x- ’ ^X- ) ~ {he S — l,he 1) , (f^xo; ^Xo ) ~ {he S — 1,/lg). 

Their conformal spins are thus s, {s — 2) and (s — 1), respectively. 

Hence, we see that we can interpret (1631) as a discrete version of the perturbed current 
conservation equation ()6.9p . with the parameters in the latter related by 

r oc ((5+ — 5^). (6.10) 
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6.3 The Ising case 

To make contact with previous work m, it is convenient to introduce the “bare” current 

J = Z®Z®---®Z®X . (6.11) 

For this operator, the linear relations (|5.1.SI) and (|5.21H arising from jgp and jen when written 
in terms of the Ising parameterisation (I2.24I[2T2^ . read respectively: 

- 0(/3.)//i(/3r.) J(F4) = 0, (6.12) 
+ QmHi{^)J{r2) - 0i(/3.)i/(/3.) J(f3) + 0i(/3,)i/(/3,) J(r4) = 0 . (6.13) 

When performing the p —?• 0 expansion using (I2.27p . the terms of order vanish in the 
combination (I6.12p — i (16.1311 . and we obtain: 

6 zi V'(ri) + 6 z 2 V’(^ 2 ) + < 5^:3 V’(^3) + i>{r4) = 

- ip[t~^'ii;{fi)+ t'ijj{f2)+ t~^'ijj{f3)+t‘ip{f4)] , (6.14) 

where t = —and we have defined 

"(/^(r) = J(r), z/i(r) = J(r). (6.15) 

In the RHS of (j6.14l) . the sum of coefficients multiplying the z/^’s is (—4zpsin0). Since the area of 
the plaquette is sin0, the relation (16.1411 is thus a discrete version of the massive Dirac equation 
with mass m = 4p r\j 

dip = —im z/i. (6.16) 

Moreover, note that at 0 = '7r/2 we have t = 1, we recover the simple form of [19| : 

6 zi z/>(ri) + 6 z 2 z/>(r2) + 5z3 ipirs) + ^^4 •0(^4) = -ip [0(n) + 0 (^ 2 ) + 0(^3) + 0(^4)] ■ (6.17) 

Finally, if we use the linear relations for the eg and ei currents instead of Bq and ei, we find the 
second part of the Dirac equations, dtp = im 0. 

7 Conclusions 

We have constructed the quasi-local operators associated to the Uq{5l2) symmetry underlying the 
chiral Potts model, for any choice of integrable Boltzmann weights. The half-currents associated 
with these operators, when dressed with suitable local phase factors, satisfy “twisted” discrete 
Cauchy-Riemann equations (15.1311 of the form 

e^^rlN^zi 0{ri) + 0 {f 2 )e-^^’^^^dz^ 0{n) + e-^’f>‘/^5z4 ©(Fi) = 0 , 

where cpr and (ps are the functions of the spectral parameters r and s along the two 

directions of the lattice. At the isotropic critical point (FZ clock model), we have exhibited the 
algebraic origin of the lattice Zw-parafermions of [18]. In the generic case N > 3, and in the 
vicinity of the FZ point, we have shown that the above equation actually encodes (a discrete 
version of) the modified current conservation relation induced by a chiral perturbation of the 
Zw-parafermionic CFT. In the Ising case, this equation also allows us to recover the discrete 
massive Dirac equation of m- 

Thus, in the framework of the chiral Potts model, we have shown that the quantum group 
symmetry can be exploited to construct off-critical discrete parafermions and to probe the nature 
of the underlying perturbed conformal field theory. 
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